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ABSTRACT 

We estimate the vertical gradient of rotational velocity for several spiral galaxies in the framework 
of a global thin-disc model, using the approximation of quasi-circular orbits. We obtain gradients 
having a broad range of values, in agreement with measurements, for galaxies with both low and 
high gradients. To model the gradient, it suffices to know the rotation curve only. We illustrate, 
using the example of galaxy NGC 4302 with particularly high gradients, that mass models of 
galactic rotation curves that assume a significant spheroidal mass component reduce the predicted 
gradient value, which may suggest that the mass distribution is dominated by a flattened disc-like 
component. We conclude that the value and behaviour of the vertical gradient in rotational 
velocity can be used to study the mass distribution in spiral galaxies. 
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1 INTRODUCTION 

Rotation curves carry important, though very limited, information 
about the kinematics of spiral galaxies. Until recently, rotation was 
determined mainly for matter moving in the direct neighbourhood 
of the galactic mid-plane. However, thanks to the constantly improv- 
ing quality of the measurements, for several galaxies it has become 
possible to dete rmine the vertical s t ructure of the rotation above the 
galactic disc (Fraternali et al. 2002; Barbieri et al. 2005; Heald at al. 
2006 j,b; iHeald et alj 12007) ; lOosterloo. Fraternali & Sancisil |2007| ; 
Levine. Heiles & B litz 2008). This enables researchers to examine the 
mass distribution in spiral galaxies more accurately. The observations 
reveal that the azimuthal component of the galactic rotation falls off 
with the vertical distance from the galactic mid-plane. The direction 
of the fall-off distinguished may indicate that the galaxy dynamics is 
governed mainly by the gravity of masses concentrated close to the 
galactic mid-plane, forming a flattened, disc-like object. As for the 
high gradient values, one can hypothesize that the flattened galaxy 
component should comprise the gross dynamical mass of the galaxy. 

Having this phenome nology in mind, we decided in our previous 
paper Jjalocha et alJl20l5T) to predict the gradient value in the Milky 
Way Galaxy. A flattened mass component was represented in our 
model by a thin disc, while the spheroidal mass components were rep- 
resented by a spherically symmetric model. We analysed how the gra- 
dient estimates depended on the relative dynamical masses ascribed to 
these components. We observed that both the value and the behaviour 
of the predicted fall-off were strongly dependent on the assumed mass 
model. We obtained correct gradient values when the gross dynami- 
cal mass was attributed to the disc component, in agreement with our 
hypothesis. We also found that similar agreement with gradient mea- 



surements was possible in the global disc model of galaxy NGC 891. 
We therefore decided in the current work to carry out similar esti- 
mates for other four galaxies (NGC 2403, 4559, 5775 and 4302) in 
which the rotation structure has been determined above the galactic 
mid-plane. 

In the current work we show that the thin-disc model, despite 
its simplicity, predicts correct values for the vertical gradient of rota- 
tional velocity, in agreement with the gradient measurements in the 
galaxies examined. In particular, this model predicts high gradient 
values that are difficult to explain by alternative gradient models (for 
example, entirely ballisti c models of disc-ha lo flow give much lower 
values than observed; see Heal d et alj J2007k 

Furthermore, the different gradient behaviour predicted by vari- 
ous mass models for the same rotation curve offers one the opportu- 
nity to test the qualitative properties of the mass distribution in spi- 
ral galaxies. High gradient values imply a flattened mass distribution 
rather than a spheroidal one. This is particularly evident when the 
gradient is measured close to the mid-plane and beyond the centre, 
further away from the galactic bulge. In this region the shape of the 
gravitational field is determined mainly by the properties of the mass 
distribution in the galactic disc and halo, while the structure of the 
bulge is not very important there, owing to the fact that the contribu- 
tion from the higher gravitational multipoles of a nearly spherically 
symmetric bulge falls off quickly with distance from the centre and 
the bulge can effectively be treated there as a point mass located in 
the centre. A qualitative argument for the suggestion that flattening 
of the mass distribution increases the gradient is provided by the fol- 
lowing reasoning. In the case of a pure thin disc, one is more likely 
to obtain high gradient values in the mid-plane vicinity than in the 
opposite case of a purely spherical mass distribution. The reason is 
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that, for the same rotation curve, the predicted gradient under the as- 
sumption of spherical symmetry is zero at the mid-plane, whereas it is 
non-zero in the disc model. Consequently, for a disc-like (i.e. mainly 
flattened) mass distribution, it is easier to produce high gradient val- 
ues close to the disc plane, whereas a feature of a nearly spherically 
symmetric mass distribution is that close to the galactic plane its con- 
tribution to the overall gradient value is low. Thus, when the mass 
model is changed so as to reduce the disc mass component at the 
cost of increasing the contribution from the spheroidal mass compo- 
nent (making the overall mass distribution more spheroidal rather than 
more flattened and keeping sufficient the amount of mass required for 
the observed rotation velocity), the predicted gradient value will have 
been decreased, especially at larger radii (where the mass distribution 
becomes in this case even more dominated by the spheroidal com- 
ponent) and closer to the mid-plane. Therefore, to ensure high gradi- 
ent values, the disc mass contribution to the mass function cannot be 
too small. An example of this interplay between various mass com- 
ponents will be discussed later in the example of galaxy NGC 4302 
with its high measured gradient. 



1.1 Vertical gradient estimates in the global disk model 

The gross dynamical mass distribution in a flattened galaxy may be 
approximately described by a thin disc, assuming axial symmetry and 
circular concentric orbits. In addition, it can be also assumed that the 
azimuthal component of the rotation above the galactic mid-plane 
can be determined by equating the radial component of the grav- 
itational force to the centrifugal force (in cylindrical coordinates). 
This is called the quasi-circular orbits approximation. Validity of 
this approximation was substantiated in the case of the Milky Way 
Galaxy by analysi ng the motions of te st bodies in the gravitational 
field of a thin disc I Jal ocha et alj|201(j|) . In spite of its simplicity, the 
model turned out to work strikingly well for the gradient description 
in the Galaxy (thi s approximation was discussed in more detail in 
Jalocha et al. 2010). We expect that this simplified description should 
be applicable in general for similar galaxies. 

In the disc model, the column mass density associated with the 
galactic mass distribution is represented by the equivalent surface 
mass de nsity of a substitute thin disc, in acc ord with the following 
formula (Jaloc ha. Bratek & Ku tschera 200~8bj) 
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Here, the integration is understood in the principal value sense (both 
integrands are divergent, however the result of integration is finite). 
This substitute surface mass density can be used for an approximate 
description of the gravitational field produced by a flattened galaxy. 
This in turn, assuming the circular orbits approximation, provides us 
with a simple tool for predicting the azimuthal velocity component 
above the galactic mid-plane and, consequently, to determine the ve- 
locity fall-off rate as a function of the vertical distance from the mid- 
plane. 

Given a surface mass density, the azimuthal component of the 
velocities of test particles in the gravitational field of a thin disc can 
be estimate d in the circular or bits approximation, using the following 
expression jjatocha et aifeoid) : 
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The model gradient value calculated as d z v v (r, z), is a function of r 
and z, whereas the gradient measurements are presented by giving a 
single number (and its error), being some average over an extended 
region. Therefore, to compare with, in a given region one could deter- 
mine a mean value of d z v v (r, z) over that region and the correspond- 
ing deviation from the mean. 

Other methods of vertical gradient estimation seem more rele- 
vant in the context of realistic gradient measurements, however. To 
mimic such measurements, one can specify in a given region an array 
of azimuthal velocities predicted by l(2j, and calculated for various 
pairs (r, z). For each specified variable z, separately, one can calcu- 
late the mean velocity by carrying out a summation over all radii in 
the region and then finding the slope of a linear fit to the mean- value 
data found at various z. We call this gradient estimation the method-I. 
Another estimation of the gradient value is obtained by finding the 
slope of a linear fit to the velocity values calculated for various z, 
separately at each fixed r, and then finding an average of the slopes 
taken over all radii in this region. This we call the method-II gra- 
dient. When gradient d z v v (r, z) is almost constant in a given region, 
both the methods should give results consistent with each other. When 
d z v v (r, z) is not constant in that region, the mean values over that re- 
gion, determined based on both methods, could be still regarded as 
describing the observed fall-off rate well, provided the dispersion of 
the slopes discussed was comparably small. 

There are other, more or less compl icated methods o f gradient 
estimation possible, which we used in djalocha et alj|20icl) . but all 
they led to similar results in the case of the Milky Way Galaxy. We 
expect that various method of gradient estimation should also give 
results consistent with each other for other spiral galaxies, therefore 
we decided here to use only the two methods described above. 



2 THE RESULTS 
2.1 NGC 2403 

We use two rot atio n curves for galaxy N GC 2403, published in 
ISofue at "3 119991 and Fraternali et al. 2002. We assume the distance 
of 3.24 Mpc, similarly as in lSofue at all 9991 . and find the integrated 
mass in disk model to be 5.22 x 10 10 Mq. Unfortunately, owing to 
the insufficiently high inclination, it was impossible to determine pre- 
cisely the vertical extent of th e gas clouds, the velo city of which was 
measured above the cold disc Fratern ali et al 1 l2002l . It was only pos- 
sible to estimate the gas-layer thickne ss to be of about 3 kpc and the 
cold-disc thickness of about 0.4 kpc dFraternali et alj |2002). There- 
fore, we assumed the region z £ (0.4, 3) kpc for the g radient es- 
timat ion. In Fig. [T] are shown th e rotation curves from (Sofue at al 
1999) and ( Fraternali et al. 2002), measured at the mid-plane. Both 
the curves are consistent with each other. Also shown is the rotation 
curve of the anomalous gas observed above the disc, somewhere in 
the interval mentioned. In between the curves are shown the rotation 
curves predicted in our model at various heights abo ve the mid-plane. 
Our results may suggest that dFraternali et alj|2002l) 's rotation curve 
of the anomalous gas, was measured even slightly above 3 kpc. 

We estimated the gradient value using the method-II in the region 
r G (1, 16.5) kpc defined by the boundary points of the anomalous 
gas measurements. We assume the radial step size of 0.5 kpc. The 
vertical interval is 2 G (0.4, 2.4) kpc, with the step size of 0.4 kpc. 
We did not use the method-I, which might have led to large errors 
owing to very extended radial region (method-I is best for the regions 
over the flat fragment of rotation curve, where the rotation above the 
disk is roughly independent of the radial variable). 

The method-II gradie nt value is — 10±4 km s _1 kpc -1 , consis- 
tent with [Fraternali (2009)'s result of about —12 km s" 1 kpc" 1 (the 
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Figure 1. Galaxy NGC2403: disc rotation curve (solid circles) and anomalous 
gas ro tation curve above the disk (open circles), both from iFraternali et"al] 
»20Q2I) . The solid lines are rotation curves predicted by our model at various 
heights above the m id-plane at 2 = .6, 1.2, 1.8, 2.4 and 3.0 kpc. The top- 
most solid line is the lSofue at a] ll999h 's rotation curve. 



Figure 3. Galaxy NGC 4559: disc rotation curve (solid circles) an d the anoma- 
lous g as rotation curve above the disk (open circles), both from lBarbieri et all 
120051) . The solid lines are rotation curves predicted by our model at various 
heights above the mid-plane at z = 0.6, 1.2, 1.8, 2.4, 3.0 and 3.6 kpc. The 
dashed line is the rotation curve predicted at z = 4 kpc. 
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Figure 2. Galaxy NGC2403. Gradient values predicted in the global disc 
model at various heights above the mid-plane. 

error was not published). Fig. [2] illustrates the gradient behaviour at 
various heights above the mid-plane. 

2.2 NGC 4559 

The vertical gradient in this galax y was observed in the region z £ 
(0.2, 4) kpc tearbierietal.ll2005l) and was determined in the same 
way as for galaxy NGC 2403. At the assumed distance of 9.7 Mpc 
to this galaxy , the i ntegrated mass corresponding in our model to 
iBarbieri et ail (2005)'s rotation curve, is 4.82 x 10 10 M Q . In Fig. 
[3] are shown the rotation curves of the cold disk of the anomalous 
gas and the rotational velocities predicted in our model at various 
heights above the mid-plane, out to z = 4 kpc, a region in which 
the rotation of gas has been observed. As previously, the exact height 
of measurements is unknown but is located somewhere in the inter- 
val z £ (0.2, 4) kpc. Our predictions are consistent with the mea- 



Figure 4. Galaxy NGC 4559: gradient values predicted in the global disk 
model at various heights above the mid-plane. 



surements; in particular, the predicted rotation at z = 4.0 kpc over- 
laps within the error bars with the rotation curve of the anomalous 
gas. This may also suggest that the gas was measured slightly above 
4.0 kpc. The method-II gradient was determined using the predicted 
rotation curves within the radial interval r £ (1.3, 19.7) kpc using 
the same steps as for the rotation measurement points, and within the 
vertical interval z £ (0.2, 3.7) kpc with the assumed step size of 
0.7 kpc. 

Method-I is not used for the same reasons as for galaxy NGC 
2403. The method-II gradient value is - 7.2 ± 2.4 km s" 1 kpc^ 1 . 
This result is lower than that given by dFratemalil f2009) of about 
-lOkms" 1 kpc" 1 (a gain, the error was not published). Since this 
result must also have an error, we are justified to say that our predic- 
tion is consistent with gradient observations. Fig.|4]shows the gradient 
behaviour at various heights above the mid-plane. 
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Figure 5. Galaxy NGC 4302: azimuthal velocity as a function of the vertical 
distance from the mid-plane. The points represent the velocity values averaged 
over the interval r £ (2.5, 6) kpc, the bars show the standard deviation in that 
interval. 



2.3 NGC 4302 



The ro tation curve of galaxy NGC 4302 was taken from lHeald et al.l 
d2007l) . The integrated mass at the assumed distance of 16.8 Mpc 
in the disk model is 5.05 x 10 10 Mq. The galaxy has a large in- 
clination; therefore the gradient measurement of the az imuthal ve- 
locity was made a nalog ously to that for NGC 891 drleald et"al1 
120071) . iHeald et alj (2007) measured the gradient in the region z £ 
(0.4, 2.4) kpc, r £ (2.5, 6) kpc and obtained the gradient value of 
about —30 km s" 1 kpc -1 (—31 ± 19.8 km s" 1 kpc" 1 in the south 
side). For the gradient calcul ation in the disk m odel, we assumed 
the same z — r region as in IHeald et al. I d2007h . with assumed z- 
step size of 0.4 kpc and r-step size of 0.5 kpc. In figure [5] the az- 
imuthal velocity averaged over r is shown as a function of the vertical 
distance. The method-I gradient is —22.7 ± 8.4 kms" kpc -1 and 
the method-II gradient is —22.3 ± 4. km s" 1 kpc" 1 . T hese values 
overlap within the errors with those for Heald et al. (2007)'s measure- 
ments. Fig. [6] shows the rotational velocity above the mid-plane pre- 
dicted in the global disc model and in a two-component model. Fig. [7] 
shows the gradient behaviour at various heights above the mid-plane. 
NGC 4302 illustrates well the influence of the spherical component 
on the gradient properties. To contrast this with the previous situation 
of purely a thin disc, consider a spherical halo with the density profile 
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comprising most of the galaxy mass, with the remaining mass at- 
tributed to an exponential thin disc. Now, the method-II gradient is 
strongly reduced to —14.8 ± 4.1 kms" 1 kpc" 1 . Fig. [8] shows the 
least-squares fit of the two-component model. Surely a mass model 
with arbitrarily assumed halo and disc-mass profile, cannot accurately 
account for the rotation curve. With a maximal halo model one can 
account for the rotation perfectly, but the corresponding method-II 
gradient value is even smaller, —11.1 ± 4.5 km s" 1 kpc" 1 . Accord- 
ingly, for mixed mass models with various relative mass profiles, the 
expected gradient value is between —11.1 km s" 1 kpc" for a more 
spherical-like mass distribution and —22.3 km s" 1 kpc" for amore 
disc-like mass distribution. To see better the influence of a heavy 
halo, let us estimate the gradient value in the mid-plane vicinity us- 
ing the rotation data at points r £ (3.5, 7)kpc in steps of 0.5kpc, 



Figure 6. The measured rotation curve of galaxy NGC 4302 (Heald et al. 
2007, solid circles), the rotational velocity above the mid-plane predicted in 
the global disk model, shown at various heights above the mid-plane in the 
interval z 6 (0.6, 3.6) kpc in steps of 0.6 kpc (solid lines), and the rota- 
tional velocity above the mid-plane predicted by the two-component model 
presented in Fig. [8] shown in the same interval at the same heights above the 
mid-plane (dashed lines). 




Figure 7. Galaxy NGC 4302: gradient values predicted in the global disc 
model at various heights above the mid-plane. 



and 2 £ (0.06, 1.06)kpc in steps of 0.2kpc. The method-II gradient 
changes only a little in the disk model —21.1 ± 6.0 km s" 1 kpc" 1 , 
but for a "disk plus halo" model it reduces substantially, to —10.4 ± 
4.2 km s" 1 kpc" 1 . The greater the radial distance, the greater is 
the halo influence and the less important the disc component. At 
small z, the gradient strongly decreases for the spherical distribu- 
tion, whereas it is still high for the disc-like distribution. It is worth 
of noting the high gradient dispersion relative to the gradient value 
— 10.4 ± 4.2 km s" 1 kpc" 1 in this region. The gradient varies more 
with r in the spherical model than in the disk model and therefore with 
increasing halo and decreasing z the gradient will be more dependent 
on r. The observations, however, do not confirm such a behaviour. 
Therefore, to see the gradient features better it would be worth re- 
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Figure 8. The least-squares fit of a two-compone nt model of galaxy NGC 4302 
considered in the text, fitted to the rotation curve jHeald et al.l2007l solid line). 
The model rotation curve is decomposed to the exponential disc (dotted line) 
and to the dark halo component (dashed line). 

peating the rotation measurements in galaxy NGC 4302 at larger radii, 
closer to the mid-plane. 

2.4 NGC 5775 

We us e the rotation curve of galaxy NGC 5775 published in llrwinl 
dl994l) . The integrated mass is 9.65 x 1O 1O M0 at the assumed dis- 
tance of 24.8 Mpc. The gradient meas urements of t he azim uthal ve- 
locity was carried out for this galaxy in lHealdatalJj2006ah in a very 
extended region covering r £ (0, 12) kpc and z £ (1.2, 3.6) kpc. 
This excludes the use of method I, since the azimuthal velocities 
at such a wide range of radii would differ too much, and the dis- 
persion of the gradient estimation might be larger than the gradient 
value. Therefore, we used only the method-II, which gives —12.0 ± 
4.3 km s -1 kpc -1 . The assumed r — z region is the same as that 
used for the measurements, and the assumed steps are A z = 0.6 kpc 
and Ar = 1.0 kpc. The measured gradient value given in lHeald at al.l 
(2006a) is 8 ± 4 km s 1 kpc 1 and is consistent within errors with 
our predicted value. Figure[9] shows the gradient behaviour at various 
heights above the disk. Fig.[l0]shows the rotational velocity above the 
mid-plane predicted in the disc model. 



3 CONCLUSIONS 

The vertical gradient in the rotational velocity of spiral galaxies can be 
satisfactorily explained in the framework of a global thin-disc model. 
The observed direction of the fall-off in rotational velocity (in the 
direction vertical to the galaxy mid-plane) and the high observed gra- 
dient values suggest that the gross mass distribution in spiral galaxies 
might be flattened, i.e. disc-like rather than spheroidal. With this hind- 
sight, assuming that the disc comprises most of the dynamical mass 
inferred from the rotation curve, we calculated the vertical gradient of 
the azimuthal velocity in the quasi-circular orbits approximation. 

We examined four spiral galaxies, obtaining gradient val- 
ues ranging from —7.2 ± 2.4 km s -1 kpc -1 to —22.7 ± 
8.4 km s -1 kpc -1 . We observed a general tendency in the gradient 
behaviour: for sufficiently large radii, the gradient values tend to be 
independent of the height above the galactic mid-plane. 

Our predictions are consistent within errors with the measured 
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Figure 9. Galaxy NGC 5775: gradient values predicted in the global disk 
model at various heights above the mid-plane. 
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Figure 10. Galaxy NGC 5775. From the top - the rotation curve from Irwin 
(19941) . rotation curves predicted in disk model at various heights above the 
mid-plane for z £ (0.6, 3.6) kpc in steps of 0.6 kpc. 

gradient values. This agreement shows that the simple model of 
a global thin disc is sufficient to account for the measured, of- 
ten high, gradient values. It is worth recalling that the global disc 
model gives the high gra dient values require d for galaxies like NGC 
4302 or NGC 891 (see IJalocha et al.ll2010l) . which are difficult to 
explain by other models such as the ballistic model, the latter of- 
ten predicting gradient values several times lower than those ob- 
served. Depending on the parametrization, the ballistic model predicts 
the gradient values of —8 km s -1 kpc -1 for NGC 4302, but more 
realistic parameter s give lower values of about —1 kms -1 kpc -1 
dHeald et al.l 120071) . much smaller than the measured ones, even 
reaching —30 kms -1 kpc -1 , whereas the averaged rotation fall- 
off we obtain for this galaxy is —22.7 ± 8.4 km s -1 kpc -1 . 
For galaxy NGC, 57 75 the ballistic mo del gives from —2.4 to 
-4.3 kms -1 kpc -1 dHeald at alj 12006a). again less than in our 
model (we obtain —12.0 ± 4.3 km s -1 kpc -1 ), but, since the mea- 
sured gradient is comparably low (—8 ± 4 km s -1 kpc -1 ), the bal- 
listic model prediction also overlaps with the measurement. However, 
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for galaxies with high gradient values our gradient modelling gives 
much better results than the ballistic model. We stress that here we 
simply compare, using the example of particular galaxies, the results 
of our model with the results of the ballistic model (with assumed pa- 
rameters similar to those in the cited papers), and we do not make any 
assessment of the ballistic model with respect to the concepts under- 
lying it. 

It is important to stress that, based only on the mass distri- 
bution precisely related to a given rotation curve, the disc model, 
apart from accounting for high gradient values, gives the correct 
low gradient when the measured gradient is low. Of course, other 
models are also capable of giving correct results for galaxies with 
low gradients. For example, in the case of galaxy NGC 5775, the 
disc model predicts —12.9 ± 2.8 km s _1 kpc -1 , the maximal halo 
model —6.5 ± 2.2 km s -1 kpc -1 , and the ballistic model gives from 
—2.4 km s -1 kpc -1 to —4.3 km s -1 kpc -1 . All of these values 
agree with the measurement — 8 ± 4. It is therefore important to have 
precise rotation measurements in the mid-plane vicinity at large radii, 
since they could discriminate between flattened and halo-dominated 
mass distributions. This is especially important for galaxy NGC 4302 
with its high gradient. 

The results of the current paper, and of our previo us paper con- 
cernin g the gradient study in the Milky Way Galaxy jjatocha et al.l 
2010), provide strong arguments that the gross mass distribution, at 
least in some spiral galaxies, might be flattened disc-like rather than 
spheroidal, unlike the suggestion of massive spheroidal dark halo 
models. Another conclusion is that the gas motion above the galactic 
discs of spiral galaxies might be governed by the gravitational poten- 
tial of the galactic disc alone rather than by non-gravitational physics 
such as matter fountains. 
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